Symplectic Integration of Constrained Hamiltonian Systems
نویسنده
چکیده
A Hamiltonian system in potential form (H(q, p) = p'M~ 'p/2 + E(q)) subject to smooth constraints on q can be viewed as a Hamiltonian system on a manifold, but numerical computations must be performed in R" . In this paper, methods which reduce "Hamiltonian differential-algebraic equations" to ODEs in Euclidean space are examined. The authors study the construction of canonical parametrizations or local charts as well as methods based on the construction of ODE systems in the space in which the constraint manifold is embedded which preserve the constraint manifold as an invariant manifold. In each case, a Hamiltonian system of ordinary differential equations is produced. The stability of the constraint-invariants and the behavior of the original Hamiltonian along solutions are investigated both numerically and analytically.
منابع مشابه
Symplectic Integration of Constrained Hamiltonian Systems by Composition Methods
Recent work reported in the literature suggests that for the long-term integration of Hamiltonian dynamical systems one should use methods that preserve the symplectic structure of the ow. In this paper we investigate the symplecticity of numerical integrators for constrained Hamiltonian systems with holonomic constraints. We will derive the following two results: (i) We show that any rst or se...
متن کاملEquivariant Constrained Symplectic Integration
We use recent results on symplectic integration of Hamiltonian systems with constraints to construct symplectic integrators on cotangent bundles of manifolds by embedding the manifold in a linear space. We also prove that these methods are equivariant under cotangent lifts of a symmetry group acting linearly on the ambient space and consequently preserve the corresponding momentum. These result...
متن کاملSymplectic Numerical Integrators in Constrained Hamiltonian Systems
Recent work reported in the literature suggests that for the long-time integration of Hamiltonian dynamical systems one should use methods that preserve the symplectic (or canonical) structure of the ow. Here we investigate the symplecticness of numerical integrators for constrained dynamics, such as occur in molecular dynamics when bond lengths are made rigid in order to overcome stepsize limi...
متن کاملSymplectic and symmetric methods for the numerical solution of some mathematical models of celestial objects
In the last years, the theory of numerical methods for system of non-stiff and stiff ordinary differential equations has reached a certain maturity. So, there are many excellent codes which are based on Runge–Kutta methods, linear multistep methods, Obreshkov methods, hybrid methods or general linear methods. Although these methods have good accuracy and desirable stability properties such as A...
متن کاملStrongly hyperbolic Hamiltonian systems in numerical relativity: Formulation and symplectic integration
We consider two strongly hyperbolic Hamiltonian formulations of general relativity and their numerical integration with a free and a partially constrained symplectic integrator. In those formulations we use hyperbolic drivers for the shift and in one case also for the densitized lapse. A system where the densitized lapse is an external field allows to enforce the momentum constraints in a holon...
متن کاملSymmetric multistep methods for constrained Hamiltonian systems
A method of choice for the long-time integration of constrained Hamiltonians systems is the Rattle algorithm. It is symmetric, symplectic, and nearly preserves the Hamiltonian, but it is only of order two and thus not efficient for high accuracy requirements. In this article we prove that certain symmetric linear multistep methods have the same qualitative behavior and can achieve an arbitraril...
متن کامل